Silica-Rich Bright Debris Disk around HD 15407A 
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ABSTRACT 

We report an intriguing debris disk towards the F3V star HD 15407A, in which an extremely 
large amount of warm fine dust (~ 10 -7 M®) is detected. The dust temperature is derived as ~ 
500-600 K and the location of the debris dust is estimated as 0.6-1.0 AU from the central star, 
a terrestrial planet region. The fractional luminosity of the debris disk is ~ 0.005, which is much 
larger than those predicted by steady-state models of the debris disk produced by planetesimal 
collisions. The mid-infrared spectrum obtained by Spitzer indicates the presence of abundant 
/mi-sized silica dust, suggesting that the dust comes from the surface layer of differentiated large 
rocky bodies and might be trapped around the star. 



circumstellar matter 



Subject headings: 
(HD 15407A) 

Introduction 



As we observed signs of col lisions of comets onto 
Jupiter in 1994 and 2009 fe.g. lHammel et alj2010h 
and c ollisions between asteroids fe.g. llshiguro et al.l 
201 if ), collisions between small bodies and/or 
planets are common in the Solar System. In 
the younger Solar System, collisions occurred 
more frequently and were one of the most im- 
portant processes in its fo rmation and evolution 
(JKenvon fc Bromlevl 120041 ) . Especially, hyperve- 
locity collisions have direct links to catastrophic 
phenomena on planets, such as craters on the 
planets, moon formation through the giant im- 
pact, and mass extinction on the Earth (e.g. 
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Bottke et al.l 120071 ). In the extrasolar system, 
products of collisions among planetesimals around 
youn g stars are observed as debris disks (jWvattl 



2008). 



HD 15407A (HIP 11696) is an F3V main- 
sequen ce star, whose dist ance from the Sun is 
55 pc (jvan Leeuwenl 120071 ). The age, the effec- 
tive temperature, and the metallicity of the star 
have been estimated as 2.1 ± 0.3 Gyr, 6350 K, 
and [Fe/H] = +0.0 8, respectively, by G eneva- 
Copenhagen Survey ( Holmberg et al.l 120091 ) . The 
parameters are confirmed by high-dispersion spec- 
tra obtained with Okayama Astrophysical Obser- 
vatory (OAO)/HIDES and the comparison of the 
position of HD 15407A in the H-R diagram to 
theoretical tracks a nd isochrones based on the 
method described in iTakedal ( 2007 ) . On the other 
hand, a recent paper suggests a much younger 
age (80^° Myr) from the lithium 671 A absorp 



tion and the motion in the Galaxy ([Melis et al 
20101 ). HD 15407A is suggested as a star in dou- 



ble system with the nearby K2V star HD 15407B 
(HIP 11692) in the SIBMAD database. The sep- 
aration between the two stars is 21'.'2 and the 
projected distance is 1170 AU at the distance of 
55 pc. 

Large excess over the expected stellar photo- 



spheric emission at 12 and 60 /im, whose color 
temperature is estimated as 750 K, has been sug- 
geste d by IRAS observations (|Oudmaiier et al.l 
1992j). The infrare d excess was confirm ed by 



the recent AKARI (jMurakami et al.l 120071 ) mid- 
i nfrared (MIR] All-S ky Survey at 9 and 18 /mi 
(jlshihara et al.l 120101 ) . suggesting the presence of 
a warm dust disk. 

In this Letter, we present new MIR photome- 
try of HD 15407A obtained with AKARI as well 
as a MIR low resolution spectrum obtained with 
Spitzer. We report the detection of abundant sil- 
ica dust together with amorphous silicate dust to- 
wards the star and discuss the nature of the debris 
disk around HD 15407A. 

2. Observations and Data Reduction 
2.1. AKARI/IRC all-sky survey 

The S9W (9 (im) and L18W (18 /im) images of 
HD 1 540 7 A were taken w ith the Infrared Camera 
(IRC; lQnakaet al.ll2007l) onboard AKARI as part 



of the All-Sky Survey observations (jlshihara et al 
2010). The star was observed on Aug 18-17, 
2006 and Feb 15-16, 2007 (UT). The fluxes at the 
two periods agree with each other within the un- 
certainty, indicating no appreciable variations in 
the flux. HD 15407A was selected as a proba- 
ble candidate of the debris disk with 18 /im ex- 
cess from the early products of the AKARI/IRC 
All-Sky Survey. The flux densities are deter- 
mined as 0.960 ± 0.031 Jy and 0.497 ± 0.021 Jy 
(AKARI-IRC-V1 J0230506+553254), in the S9W 
and L18W bands, respe ctively, in the AKAR I/IRC 
Point Source Catalog ( Ishihara et al.l 120101 ). Al- 
though the AKARI flux densities of HD 15407A 
might be contaminated by HD 15407B due to 
the large beam sizes (5'.'5 and 5'.'7 for the S9W 
and L18W bands, respectively), the fractional flux 
densities of the contamination are estimated to be 
< 20% in both bands and thus the detection of 
the excess is secure. 



2.2. Spitzer/TRS observations 

Follow-up spectroscopic observations of HD 15407A 
were made with t he Infrared Spectrograph (IRS; 
iHouck et al1l2004l ) onboard Spitzer on Oct 9, 2008 
(AOR ID# 26122496). All the four low-resolution 
modules, Short-Low 2 (5.2-7.7 /im) and Short- 
Low 1 (7.4-14.5 /im), Long-Low 2 (14.0-21.3 /im) 



and Long-Low 1 (19.5-38.0 /im), were used to ob- 
tain a full 5-35 /im low-resolution (A/AA ~ 100) 
spectrum. We used the pipeline-processed (S18.1) 
Basic Calibrated Data products for the target and 
analyzed the sky-subtracted extracted ID spectral 
data since the target is a point source in a rela- 
tively empty field. The wavelength calibration is 
as good as 0.1 /im along the dispersal direction. 
For the Long-Low 1 spectra, we use the data only 
for A < 35 /im because the noise becomes large at 
A > 35 /im, whose spectral range is not critical for 
the present analysis. We adopt the flux densities 
of the pipeline-processed products, for which the 
absolute accuracy in the Short-Low and Long-Low 
spectra is better than 10%. The Spitzer spectrum 
does not include the contribution of HD 15407B 
since HD 15407A and HD 15407B are clearly sep- 
arated in the peak- up image of the IRS. 

3. Results 

3.1. Spectral Energy Distribution 

The observed flux densities of HD 1540 7 A in 
various bands are shown in Table [T] The photo- 
spheric flux densities of H D 15407A are estimated 
from the Kurucz model (KuruczJ|l992) with the 
effective temperature of 6500 K and the surface 
gravity of log.g = +4.0 fitted to the 2MASS K s - 
band photometry of the star. The resultant pho- 
tospheric spectrum fitted with the Xs-band is in 
agreement with the 2MASS J- and -ff-band data. 
Therefore the extinction to HD 15407A is negligi- 
ble in the infrared. The extinction in the optical 
region is also negligible since the observed B — V 
color of the star is consistent with the intrinsic 
color of F2-5V stars (B - V = 0.35 for F2, and 
0.44 for F5 stars). 

The obtained spectral energy distribution 
(SED) of the star compiled with the Spitzer/IRS 
spectrum and the AKARI and IRA S photometry 
is shown in Figure [1] together with the expected 
photospheric emission. Significant excess emission 
at wavelengths longer than 5 /im is clearly seen 
and the detected flux densities at 9 and 18 /xm 
are 5 and 10 times larger than the photosphere, 
unambiguously indicating the presence of a warm 
and bright debris disk. 



3.2. Dust Features 

To examine the excess emission in detail, we 
subtract the estimated photospheric emission from 
the observed IRS spectrum (Figure [5]). The 
photosphere-subtracted spectrum shows two sig- 
nificant emission features centered at around 9 and 
20 /mi. A weak feature at 16 /mi is also seen. A 
ubiquitous dust species family, silicate, which has 
spectral features around ~ 10 /mi and ~ 20 /mr, 
can be considered as a main carrier of the observed 
features at the first sight. However, sub-/j,m-sized 
amorphous olivine and pyroxene commonly seen in 
celestial objects have features around 9.3-9.7 /xm 
and 18.0 /im, which cannot reproduce the ob- 
served spectrum completely. As an alternative 
carrier of the observed features around around 9 
and 20 /im, we consider amorphous silica (silicon 
dioxide), which also has Si-0 and O-Si-0 modes 
and shows broad features at similar wavelengths 
to amorphous silicate. 

3.3. Spectral Fitting 

To estimate the dust temperature and the 
species, we perform a fit with an SED model that 
consists of blackbody emission, silicate and silica 
emission. The model flux density is given by 

^exc^M = ^BB-Bi/(A,T B b) 

i ^silicatc^silicatc-^'^ (,A, 1 silicate J 
i ^silica^silica-t'^ \" 7 -L silica J > 

where J^bb is the solid angle of a blackbody rep- 
resenting the continuum emission from large dust 
(zL 10 /mi in size), B„(A,T) is the Planck func- 
tion with the temperatures T, ft (A) is the mass 
absorption coefficients of the dust, and a's are the 
scaling factors proportional to the dust mass. For 
silicate, we choose one from eight mass absorp- 
tion coefficients (amorphous olivine or pyroxene 
with the size of 0.1, 1.0, 1.5, or 2.0 /xm) com- 
puted from the o ptical constants o f MgFeSiC>4 and 
Mgo.sFeo.sSiOs (Dorschne r et al.1 1995[) based the 
Mie theory (Bohr en &: Hufmannlll983f ). We take 
the specific mass densities of both olivine and 
pyroxene as 3.3 g cm" 3 . For silica, we choose 
one from fuse d quartz (IKoike et alJll989l ) or an- 
nealed silica ( Fabian et al.l 12000 ). Fused quartz 
is an amorphous (glassy) polymorph of silica and 
shows broad features in the MIR. Annealed sil- 
ica may contain cristobalite and tridymite, poly- 



morphs of crystalline silica, and shows relatively 
sharp features. We take the mass absorption co- 
efficients of annealed silica and compute that of 
fused quartz for the shape of con tinuous ellipsoidal 
distribution (JFabian et al.l 120011 ) fro m the optical 



const ants measured in laboratories ([Koike et al 
19891 ). We search for parameters minimizing the 



reduced xt using the data between 5 and 35 /mi 
with respect to each of Obb and a's. 

At first, we find that a combination of black- 
body, 1.5 /mi-sized amorphous pyroxene and fused 
quartz provides the best fit (Figure[2]) . The second 
best combination is blackbody, 1.5 /mi-sized amor- 
phous pyroxene, and annealed silica. The best- 
fit model spectra and parameters of these combi- 
nations are shown in Figure [2] and Table [2] re- 
spectively. The models with the two combina- 
tions of dust components seem fairly good, sug- 
gesting that it is a robust conclusion that sil- 
ica dust is one of the major dust species of de- 
bris dust together with amorphous pyroxene dust 
around HD 15407A. The model spectrum with an- 
nealed silica explains the observed small feature 
at A = 16 /zm while it produces a sharp extra 
emission at A = 12.6 /mi, which is not observed. 
Since fused quartz and annealed silica have similar 
but slightly different band features, we perform a 
tuned fit assuming that the silica component is a 
mixture of fused quartz and annealed silica, and 
that the silicate is 1.5 /zm-sized amorphous py- 
roxene with T s iii ca t = Tallica. The best-fit model 
spectrum and parameters in this dust combination 
are shown in Figure [2] and Table [2j respectively. 
Since this model provides the least xt among 
the models considered here, we adopt this model 
(T BB = 505 K and T si i icate = T siVlca = 615 K) 
in the following discussion. The masses of amor- 
phous pyroxene and silica dust around HD 15407A 
are derived as 4.2 x 10 17 and 2.5 x 10 17 kg, respec- 
tively, assuming that the dust emission is optically 
thin. The derived total mass of the fine dust in the 
warm debris disk is 6.8 x 10 17 kg - 10~ 7 M e . This 
value does not include the mass of the blackbody 
dust, which is supposed to come from large-sized 
(> 10 /mi in size) grains and rubbles. 



4. Discussion 

4.1. Radial Distribution of Dust 

We estimate the distances from the central star 
of the blackbody and the 1.5 /mi-sized amorphous 
pyroxene dust from the dust temperatures and 
emissivities assuming that the grains are in ra- 
diative equilibrium around a F3V star with the 
luminosity of 3.92^0 . The derived distance for 
the blackbody and the pyroxene dust is 0.6 and 
1.0 AU, respectively, which corresponds to a ter- 
restrial (rocky) planet region. The silica dust must 
contact with other dust species to keep its tem- 
perature high since it is transparent in the visible 
wavelength and the stellar radiation solely can- 
not heat silica grains to the observed temperature. 
Thus, although the distance of silica from the star 
in radiative equilibrium is < 0.1 AU, the silica dust 
is considered to be located at a region similar to 
those of pyroxene and blackbody dust. Only upper 
limits are available for the far-infrared (FIR) flux 
densities from IRAS observations and the exten- 
sion of the debris disk cannot be estimated. High 
spatial resolution observations and FIR-radio ob- 
servations are required for further examination of 
the radial distribution of dust around the star. 

4.2. Fractional Luminosity 

The fractional luminosity (/ obs = L dust /L stllI ) 
of the warm debris disk is estimated as ~ 0.005 
from the fitting of the SED. Figure 3 plots the 
maximum fractional luminosity / max predicted by 
a simple model of the steady-state evolution of 
debri s disks produced by collisions ( Wvatt et al.l 
20071 ) with the stellar parameters of HD 15407A 
and the derived dust disk radii (i?dust = 0.6 and 
1.0 AU for blackbody and pyroxene dust, respec- 
tively). The ratio of /obs//max of the debris disk 
around HD 15407A is larger than 10 5 if we as- 
sume 2.1 Gyr as the age o f HP 1540 7 A as sug- 
gested by iHolmberg et al.l (|2009l ) . The observed 
fractional luminosity of HD 1540 7 A system cannot 
be accounted for by a simple steady-state model. 
Even whe n we assume the a ge of 80^20 Myr sug- 



gested by iMelis et al.l (|2010l ). / bs//max is about 



10 and significant enhancement of / Q bs compared 



to f„ 



secure. To date several warm de- 



be responsible; e.g. IWvatt et al.l ([20071 ) estimated 
W/max = 10 3 -10 4 for BD +20 307, HD 72905, 
T) Corvi, and HD 69830, which possess warm dust 
at 0.2-2 AU. HD 15407A is one of the highest- 
/obs//max debris disks among these and seems a 
non-typical warm debris disk. The disk is not a 
protoplanetary one since H/3, H7, and US lines 
are seen in absorption, not in emission, in the 
OAO/HIDES optical spectrum of the star, which 
suggests the absence of gas accretion around the 
star. 

HD 15407 is a possible binary system and 
HD 15407B might disturb circumstellar material 
around HD 15407A dynamically. At present there 
is no theoretical model of debris disk evolution 
for a binary system in terms of fractional lumi- 
nosity. However, the study of the stability zone 
in a binary system by iHolman fc Wiegerd ( 19991 ) 
suggests that material at < 400 and < 50 AU 
from HD 15407A is dynamically stable for cases 
where the eccentricity of the system is 0.0 and 0.8, 
respectively. 

The radiation pressure should blow out /un- 
sized small grains, which show prominent feature 
in the MIR, in the vicinity of a F3 star in a 
short (< 100 years) time scale. It cannot be con- 
cluded only from available data whether or not 
HD 15407A shows a time variation of its fea- 
ture emission. Time variation in the continuum 
flux level of the excess is not seen in the ob- 
servations after IRAS, which can be accounted 
for since grains responsible for the continuum ex- 
cess are large (> 10 /xm) and thus have long 
(> 1000 years) lifetime. Monitoring observa- 
tions of MIR spectrum would be important to 
examine the actual lifetime of small dust around 
HD 15407A. If dust features towards HD 15047A 
do not change over decades, mechanisms that re- 
tain the fine grains near the star should be playing 
a role in HD 15407A, such as shepherding of dust 
aroun d a planet as in Saturn (iSpahn fc Sponholz 
19891) and Uranus (|Murrav fc Thompsonlll990l ) or 



bris disks are known, whose fractional luminosi- 
ties are much larger than the steady-state model 
and in which transient events are suggested to 



dust capt ure in a resonance as in the circumsolar 
dust ring ( Reach et al.lll995l ). which are in effect 
in the Solar System. Alternatively the fine dust 
may be produced violently and continuously over 
decades around HD 15407A . It should be noted 
that the IWvatt et al.l ( 20071 ) model does not as- 
sume the presence of planet and trapping dust in 
a resonance in a planetary system. There is no 



evidence of a planet orbiting around HD 15407A 
at present. But if trapping of dust takes place ef- 
ficiently, the predicted fractional luminosity from 
the model may be different than in debris disks 
without planets. 

4.3. Silica in Debris Disks 

The mineralogical characteristics of the de- 
tected dust provide a clue for the origin of the 
debris dust. In the Solar System, silica is one of 
the most abundant minerals in the Earth's crust. 
The presence of amorphous silica is also sug- 
gested from the analysis of the Comet 81P/Wild 2 
dust sample retrieved by the STARDUST mission 
( Mikouchi et al.l 120071) . The presence of annealed 
silica dust in the protoplanetary disk around sev- 
eral T Tauri stars (TTSs) has been indicated i n 
their Spitzer MIR spectra \ Sargent et al.l 120091 ) . 
However, silica is not a major constituent of the 
interstellar medium since the expected 9 /im fea- 
ture is not seen in the interstellar medium towards 



the Galactic Center ([Kemper et al.l 12004). 

The presence of silica-like dust has also been 
suggested in debris disks. An almost equiva- 
lent amount of /im-sized amorphous pyroxene and 
fused quartz are detected towards HD 15407A in 
addition to the blackbody component. Ground- 
based observations of 8-13 /xm spectroscopy in- 
dicate the presence of the 9 /xm-peaked broad 
feature which is attributable to amorphous silica 
around the 100 M yr-old GO field star HD 23514 
(jRhee et al.l 120081 ). Spitzer/IRS spectroscopy re- 
veals the presence of abundant amorphous silica 
dust around the A0V st ar HD 172555 who se age is 
estimated as ~ 12 Myr ( Lisse et al.ll2009l ). Detec- 
tions of silica among protoplanetary disks around 
TTSs, debris disks, and the Solar System suggest 
that silica might be ubiquitously present through 
the planet formation processes. Debris disk stars 
with silica dust share a common habit that they 
show excesses emission in the near-infrared (A < 
5 /mi) , suggesting the presence of very warm dust 
in the vicinity of the stars. The presence of very 
warm dust might be important for silica grains, 
which are almost transparent in the optical wave- 
length and thus needed to be heated by contact 
with other dust species thermally and emit in the 
MIR. 

It is widely accepted that debris dust is pro- 
duced from collisions of planetesimals. There- 



fore, if a large amount of silica-rich planetesimals 
around HD 15407A collide with each other ac- 
tively, the presence of abundant silica dust around 
HD 15407A could be accounted for. Although 
the origin of silica-rich planetesimals is still an 
open question, a breakup of a large differentiated 
rocky body may be a likely possibility. A similar 
mechanism has been proposed as an origin of the 
enstatite-rich bright d ebris disk around HD 165014 
(|Fuiiwaraet alj|2010l) . 

As a probable source of silica dust around 
HD 172555, in which abundant SiO gas is also 
detected, giant hyper velocity (> 10 km s _1 ) im- 
pact bet ween large r ocky planetesimals is sug- 
gested bv lLisse et al.l ( 20091 ) . It may be similar to 
the one that formed the Moon (giant impact) or 
the one that stripped the surface crustal material 
off of Mercury's surface. The comparable silica- 
to-silicate ratio in /im-sized small dust around 
HD 15407A suggests a similar origin of the sil- 
ica dust with that of HD 172555. This also seems 
harmonic with the large fractional luminosity of 
HD 15407A's debris, which may be connected with 
a transient event. However, the presence of SiO 
gas around HD 15407A is not confirmed in its 
Spitzer/IRS spectrum. Assuming that SiO gas 
is really absent arou nd HD 15407A, a substantia l 
time (~ 10 3 -10 4 vr; IPahlevan k Stevenson] 120071 ) 
should have elapsed after the impact and SiO gas 
should have disappeared due to the process of re- 
condensation even if a hypervelocity impact might 
have occurred and have produced silica-rich dust 
around HD 15407A. Taking account of the short 
lifetime of ^m-sized dust, dust trapping and/or 
continuous dust production by collisional grinding 
of the large rubbles are required after the possi- 
ble impact. Further theoretical studies in addi- 
tion to observations with high spatial resolution 
and wide wavelength coverage will provide hints 
on the possible linkage between the large infrared 
excess, the dynamical evolution, and the possible 
transient event of the intriguing debris disk around 
HD 15407A. 
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Table 1: Infrared Photometry of HD 15407A. 
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0.96 ±0.03 


AKARI/IRC 


0.22 


25 
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1.05 ±0.06 
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0.12 


16 
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0.50 ±0.02 
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0.06 


22 
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0.43 ±0.03 


IRAS 


0.03 
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GO 


<0.40 


IRAS 


0.005 
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< 1.14 


IRAS 
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a From Kurucz model to fitted to 2MASS JHK B -bands data. 
b X = (Observed — Kurucz)/noisc. 



Table 2: Best-fit Parameters of Spectral Fitting 3 





Blackbody 


Silicate 




Silica 




Silica Component 


Tbb 
(K) 


(10~ 16 Str) 


-* silicate 

(K) 


Mass 
(10 17 kg) 


J- silica 

(K) 


Mass 
(10 17 kg) 


xi 


Fused+ Annealed 


505 


2.3 


615 b 


4.3 


615 b 


1.3 + 1.2 


5.4 


Fused Quartz 
Annealed Silica 


505 
515 


2.3 

2.4 


725 
815 


3.0 
2.5 


605 
600 


2.2 

2.6 


6.0 
6.4 



"The number of free parameters is 6 for all of the listed models. 
6 Assumed T.iij-,,4.. = T.\\\ r! ,. 
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Fig. 1.— SED of HD 15407A. The blue 
filled squares, green open circles, and purple 
crosses indicate the photometric data obtained 
with AKARI/IRC, IRAS, and 2MASS, respec- 
tively. The horizontal bars of the photometric 
data show the band widths. The red and gray lines 
indicate the Spitzer/IRS spectru m and the pho - 
tospheric contribution of F3 star ( Kurucall992 K 
respectively. 
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Fig. 2. — Debris disk spectrum of HD 15407A (red 
crosses) and the best-fit model spectrum (solid 
black line) consisting of blackbody dust (gray), 
1.5 /im-sized amorphous pyroxene (green), fused 
quartz (blue), and annealed silica (purple). The 
debris disk spectrum is derived by subtraction of 
the estimated photospheric contribution from the 
observed IRS spectrum. The residual spectra sub- 
tracted by the best-fit model is shown below the 
fit result. The model spectrum with fused quarts 
or annealed silica only as silica dust is also shown 
shown in the main panel as dashed and dotted 
black lines, respectively. 
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Fig. 3. — Fractional luminosity of the debris 
disk around HD 15407A. The theoretical maxi- 
mum fractional luminosities of planetesimal belts 
at i?dust = 0.6 (dashed line) and 1.0 AU (solid 
line) around an F3V star are calculated by a model 
of the steady-state e volution of a debris disk pro- 
duced by collisions (|Wvatt et al.1 120071 ) with the 
fixed model parameters (belt width: dr/r = 0.5; 
planetesimal strength: Q* D = 200 J/kg; planetesi- 
mal eccentricity: e = 0.05; diameter of the largest 
planetesimal in cascade: D c = 2000 km). The ob- 
served fractional luminosity towards HD 15407A 
(filled squares) , which is shown by a ssuming the 
age of 2.1 Gyr fflolmberg et alJ l2009h and 80 Myr 
(JMelis et alJl2010h . is 10 4 -10 5 times larger than 
that expected from the model. 



